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The  strain  and  electric  fields  present  in  free-standing  AIGaN/GaN  slabs  are  examined  theoretically 
within  the  framework  of  fully  coupled  continuum  elastic  and  dielectric  models,  Simultaneous 
solutions  for  the  electric  field  and  strain  components  are  obtained  by  minimizing  the  electric 
enthalpy.  We  apply  constraints  appropriate  to  pseudomorpliic  semiconductor  epitaxial  layers  and 
obtain  closed-form  analytic  expressions  that  take  into  account  the  wnrtzite  crystal  anisotropy.  It  is 
shown  that  in  the  absence  of  free  charges*  the  calculated  strain  and  electric  fields  are  substantially 
da  Herein  from  ihuse  obtained  using  the  standard  model  without  electromechanical  coupling.  It  is 
also  shown,  however,  that  when  a  two-dimensional  electron  gas  is  present  at  the  AIGaN/GaN 
interface,  a  condition  that  is  the  basis  for  heterojunction  field-effect  transistors,  the 
electromechanical  coupling  is  screened  and  the  decoupled  model  is  once  again  a  good 
approximation.  Specific  cases  of  these  calculations  corresponding  to  transistor  and  superlattice 
structures  are  discussed.  ©  2003  American  Institute  of  Physics.  [DGl:  1 0.1063/ L  1.620378] 


I.  INTRODUCTION 

In  previous  analyses  of  the  strain  in  epitaxial  layers  of 
AIGaN/GaN,  the  electrical  and  mechanical  properties  of  the 
crystal  have  been  treated  as  though  they  are  decoupled.  Lin¬ 
ear  elastic  theory1  is  assumed  to  hold,  and  Hooke's  law  is 
invoked  to  obtain  the  relation  between  the  in-plane  and  axial 
component;  of  the  strain  tensor.  The  standard  model  de¬ 
couples  the  strain  tensor  from  the  electric  field  and  enables 
the  separation  of  the  electric  field  and  electronic  eigenstate 
Calculations  from  calculations  of  the  strain  field.  Historically, 
there  is  a  solid  tradition  of  using  this  separability  when 
studying  hctcrostracturcs  of  GaAs  and  associated  alloys 
where  the  approximation  that  the  strain  field  is  negligibly 
affected  by  the  electric  field  is  valid.1 

Generally,  it  is  known  from  thermodynamics2  that  the 
electrical  and  mechanical  properties  of  piezoelectric  materi¬ 
als  are  coupled  and  a  simultaneous  treatment  is  called  for 
when  the  electromechanical  coupling  is  significant-  This  is 
Ihe  case  when  the  piezoelectric  response  is  large,  as  it  is  in 
AlGaN,  especially  for  high  Al  fractions.  For  example,  large 
corrections  of  the  electrostatic  and  elastic  properties  have 
been  predicted  for  AIGaN/GaN  transistor  structures1  and  tor 
free-standing  plates  of  AIN4  when  a  fully  coupled  model  is 
used  instead  of  a  standard  (uncoupled.)  one. 

The  electromechanical  coupling  is  the  basis  for  surface 
acoustic  wave  (SAW)  devices  using  AIN  and  GaN  thin 
films.5-*1  A  quantitative  measure  of  the  interaction  between 
the  electrical  and  mechanical  properties  of  the  material  is  the 
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electromechanical  coupling  coefficient*10  an  important  pa¬ 
ra  ineter  in  the  design  of  SAW  filters.  This  quantity  is  related 
to  the  interaction  between  acoustic  and  electromagnetic 
waves  in  piezoelectric  materials11  and  is  derived  from  the 
thermodynamic  equations  of  state  for  piezoelectric 
materials.10  The  electromechanical  coupling  described  in  the 
present  work  addresses  the  interaction  between  the  static 
electric  and  mechanical  strain  fields.  To  date,  this  interaction 
has  been  ignored  in  file  modeling  and  design  of  AIGaN/GaN 
heterojunction  field-effect  transistors  (HFETs).  The  electro¬ 
mechanical  coupling  described  herein  is  not  directly  related 
to  the  electromechanical  coupling  coefficient  in  SAW  de¬ 
vices,  but  both  originate  from  the  same  thermodynamic 
equations  of  state,10 

In  this  article,  we  calculate  the  strain  field  in  a  free¬ 
standing  AIGaN/GaN  slab  using  a  fully  coupled  model-  A 
free-standing  bi layer  slab  is  chosen  as  the  mode!  structure  on 
which  to  develop  the  theory  as  it  can  serve  as  a  building 
block  for  other,  more  complicated,  structures.  For  instance,  it 
can  be  used  as  a  building  block  for  superlattices  and  quan¬ 
tum  wells.  In  addition,  a  HFET  is  just  a  special  case  of  a 
hi  layer  slab  in  the  limit  that  the  GaN  layer  is  much  thicker 
than  the  AlGaN  layer.  It  will  be  seen  that  the  fully  coupled 
strain  and  electric  fields  for  an  HFET  can  be  obtained  in  this 
asymptotic  limit. 

Although  the  present  model  is  still  within  the  framework 
of  continuum  mechanics,  it  goes  beyond  the  standaad  con¬ 
tinuum  elastic  theory  used  in  typical  strain  calculations  in 
semiconductor  materials.  We  obtain  the  total  electric  en¬ 
thalpy  for  the  bilayer  slab  and  minimize  it  subject  to  certain 
constraints  to  find  both  the  strain  and  electric  fields.  Two 
constraints  are  used.  One  is  that  the  two  layers  must  share  a 
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common  r-plane  lattice  constant  after  strain.  This  lattice  con¬ 
stant.  is  unknown  at  the  outset  and  is  deduced  only  after 
minimization.  The  other  constraint  is  set  up  by  the  electro¬ 
static  potential,  and  lienee  the  electric  field,  being  forced  to 
satisfy  the  Poisson  equation  subject  to  the  boundary  condi¬ 
tions  at  the  surface  and  the  common  interface-  Both  the  spon¬ 
taneous  and  piezoelectric  polarizations  are  included.  The  ef¬ 
fect  on  the  strain  and  electric  fields  of  a  two-dimensional 
electron  gas  (2DEG)  at  the  AIGaN/GaN  interface,  a  situation 
that  is  essential  for  channel  conduction  in  the  HFEX  is  also 
investigated-  The  model  products  expressions  for  the  strain 
tensor  and  the  electric  field  in  the  two  layers  in  closed  form. 

This  article  is  organized  as  follows:  In  Sec.  II,  the  elec¬ 
tric  enthalpy  for  the  AlGaN/GaN  slab  is  minimized  subject  to 
the  constraints  derived  from  the  pseudomorphie  strain  con¬ 
dition  and  from  the  Poisson  equation.  The  strain  tensor  and 
the  electric  held  are  derived  in  closed  form.  In  Sec.  Ill,  the 
calculated  results  from  the  fully  coupled  model  are  con¬ 
trasted  with  those  from  the  uncoupled  model-  The  effect  on 
the  strain  and  electric  fields  of  screening  from  the  2 DEG  is 
discussed.  The  results  are  summarized  in  Sec.  IV. 

II.  MODEL  DESCRIPTION 

In  the  standard  strain  theory  for  semiconductor  materi¬ 
als,  the  formal  way  to  calculate  the  strain  held  in  a  general¬ 
ized  strain  problem  is  ro  minimize  the  Helmholtz  free 
energy 12,1  J 

F=  \Ctjuy  1)7*1’  (IJ 

in  which  C ^  is  the  fourth- ranked  elastic  stiffness  tensor,  yi{ 
is  the  strain  tensor,  and  the  indices  ri  /,  k,  and  f  run  over  the 
Cartesian  coordinates  x,  y>  and  z.  Summation  over  repealed 
indices  is  implied  throughout-  The  minimization  will,  of 
course,  be  subject  to  certain  constraints  brought  about  by  the 
boundary  conditions  at  the  surfaces  and  between  adjacent 
materials,  as  well  as  any  external  force s. 

In  a  wide  range  of  semiconductor  problems,  however, 
one  is  concerned  with  calculating  the  strain  in  epitaxial  lay¬ 
ers  grown,  at  least  in  principle,  pseudomorphieally  on  a  thick 
buffer  layer.  Typically,  the  thick  layer  belongs  to,  or  at  least 
closely  resembles,  the  same  crystallographic  point  group  as 
the  epilayers.  In  this  case,  a  minimization  of  Eq.  (1)  is  un¬ 
necessary.  since  we  can  take  full  advantage  of  the  boundary 
condition  for  a  free  surface  which  requires  that  the  compo¬ 
nents  of  flie  stress  tensor  parallel  to  an  outward  normal  to  the 
surface  be  zero.  For  a  surface  normal  to  the  z  axis,  this 
means  that  ah=Ot  where  is  die  stress  tensor.  In  addition h 
since  the  lateral  extent  of  the  layers  is  far  greater  than  their 
thickness,  a  one-dimensional  (ID)  approximation  can  be 
used,  wherein  y^-0  for  j.  This  condition  is  equivalent  to 
ignoring  the  effects  of  bowing  on  die  local  lattice  displace¬ 
ment-  Under  these  conditions,  the  result 

"z-= o  s  cxxa{  yfJ) + (2) 

is  obtained  for  wurtzite  epilayers  oriented  in  the  [0001]  di¬ 
rection.  Here  the  only  unknown  is  y.5T  hecause  a  standard 
approach  is  to  assume  that  the  buffer  layer  is  unstrained  and 
that  the  epilaycr  assumes  the  r-plane  lattice  constant  of  (he 


Fin.  t.  Cross  section  of  &  model  cation-fa  wd  AlCraK/Oah'  bilayer  slab.  The 
top  panel  shows  the  direction  of  Lhe  piezoelectric  and  spontaneous  Fs 
pojflii7itk)n  vectors  in  relation  10  U»  z  axis.  Both  layers  in  the  sLab  are 
under  strain,  and  the  directions  of  she  Four  strain  components  are  indicated. 
tB  and  tb  are  the  thicknesses  of  the  AlGaN  and  G aN  layers,  respectively  The 
bottom  panel  shows  the  fixed  space  chaises  sir  the  extremities  of  each  layer 
resulting  From  potanzahon.  AiSO  shown  in  the  boltom  panel  Ls  She  potential 
taken  to  be  the  same  on  both  surfaces,  L.c-,  the  applied  bins  is  taken  to  be 
zero. 


buffer,  effectively  fixing  yXl  and  yp.,  Clearly,  when  y^ 
=  Yyy=(a-a0)/ae>  Eq.  (2)  yields  the  usual  result  for  the 
Poisson  effect  in  the  case  of  biaxial  strain  in  an  epi layer 
which  is  y„/r^=-2C#J(-JC-^, 

The  situation  for  the  free-standing  bilayer  slab  depicted 
in  Fig.  1  in  which  the  two  layers  are  of  comparable  thick¬ 
nesses  is  more  complicated,  since  none  of  the  strain  compo¬ 
nents  is  known  beforehand.  To  obtain  the  strain,  as  least 
within  the  standard  model,  the  total  mechanical  strain  energy 
for  the  slab  must  be  minimized  with  respect  to  subject  to 
the  constraints  along  the  interface.  There  is  a  further  compli¬ 
cation  in  piezoelectric  materials  such  as  AIGaN  and  GaJN 
wherein  the  electrical  and  mechanical  properties  are  coupled 
through  the  piezoelectric  coefficient  tensor.  Equation  (1)  is 
no  longer  a  suitable  energy  functional  for  calculating  the 
strain  held,  since  it  does  not  include  the  electromechanical 
coupling.  Instead,  we  begin  the  derivation  with  the  electric 
enthalpy  H  given  hyH 

H=U- E*JX  (3) 

where  E  and  D  are  the  electric  field  and  electric  displace¬ 
ment,  respectively,  and  U  is  the  total  internal  energy  given 
by 

U  =  j  CjjuYi]  i  je  tjEsEj  t  (4) 

in  which  Rf}  is  the  tensor  form  of  die  electric  permittivity. 
The  electric  displacement  field  components  are  given  by 

Di~eijl;Yjk+eijEj+rf  ■  (5) 

s 

in  which  is  the  piezoelectric  coefficient  tensor  and  n 
the  spontaneous  polarization.15  The  first  term  in  Eq,  (5)  is 
recognized  as  the  piezoelectric  polarization.  The  spontaneous 
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polarization  is  in  the  direction  M— *Ga  along  the  oaxis  bond. 
After  substitution  into  Eq.  (3.K  the  enthalpy  in  its  final  form 
is 

^  ¥ij  ¥tt  ~  eijkE\ 7j*  ~  2  E \  j^i ~  E i^i  1  ( ti) 

This  expression  includes  the  full  elec tromechan Seal  coupling 
as  well  as  die  spontaneous  polarization.  Differentiating  the 
enthalpy  with  respect  to  the  strain  tensor  gives  the  fully 
coupled  equation  of  state 


&ij  gtijEk 

for  piezoelectric  materials.  To  obtain  the  strain  and  electric 
fields  within  a  fully  coupled  model,  Eq .  (6)  must  he  mini¬ 
mized  with  respect  to  both  yi{  and  Er-  subject  it)  the  con¬ 
straints  to  be  discussed  shortly  for  the  hilayer  slab  of  Fig.  L 
Using  the  Voigt  notation11  and  expanding  Eq.  (6),  the 
enthalpy  for  a  wurtzite  crystal  with  the.  [OGGI]  axis  as  the 
principal  axis  can  be  written  as 


H  =  TCI  |(  7^>+  ?  C  + C  isr„r.,jr +C  l3yzz(  r„+  yyf)Jr2CM(  y$t  +  yj() +{Cn  -  C  ia)  -  £.|>ii(  yxx  +  yyy) 

+  'J3T«]-*i siExyxz+E,yyt)-EtPs-  (8) 


r 


where  e  is  the  electric  permittivity  assumed  here,  for  sim¬ 
plicity,  to  he  a  scalar.  The  enthalpy  for  the  slab  is  given  by 


Hmb= 


taH*+lhHh 

G+4 


(9) 


where  t(i  and  ^  are  the  thicknesses  of  the  two  layers  and  jLaM 
denotes  the  AlGaN  layer  and  “b"  the  GaN  layer.  It  is  noted 
that  the  elastic  and  piezoelectric  coefficients  and  the  electric 
permittivity  can  be  substantially  different  in  the  two  layers. 
These  variations  are  taken  into  account  in  Eq.  (9),  Because  of 
the  rotational  symmetry  of  the  slab,  y>lv=  yxx  and  yy-  =  ysz 
in  each  layer.  The  shear  terms  in  the  strain  tensor,  y^  for  t 
turn  out  to  be  zero  because  of  the  assumption  that  the 
strain  is  piecewise  homogeneous,  A  more  realistic  assump¬ 
tion  would  be  to  expect  the  strain  to  diminish  away  from  the 
interface,  resulting  in  a  nonzero  yx~  that  would  vary  with 
position  along  the  c  axis.  This  situation  would  be  manifested 
in  a  bowing  of  the  slab.  Such  inhomogeneous  strains,  how¬ 
ever,  would  require  a  three-dimensional  OD)  numerical  cal¬ 
culation  that  is  beyond  the  scope  of  the  present  work. 

One  of  the  minimization  constraints  on  Eq.  (9)  is  that  the 
two  layers  share  the  same  c -plane  lattice  constant  which  is 
unknown  at  the  outset.  This  condition  is  expressed  as 

{J+Ti>a-(l  +  1&>a*,  (10) 


w  here  and  are  the  unstrained  c-plane  lattice  of  the 
AlGaN  and  GaN  layers,  respectively.  Equation  (10)  pre¬ 
sumes  ideal  growth  conditions.  Partial  relaxation  due  to  dis¬ 
locations  will  be  the  subject  of  future  investigations. 

The  other  constraint  is  the  relationship  between  the  elec¬ 
tric  and  strain  fields.  This  relationship  is  established  in  the 
present  work  by  solving  the  ID  Poisson  equation  given  by 


=0  and  z  =  (ffl+rir)  and  also  to  the  continuity  of  ^and  the 
electric  displacement  across  the  interface.  The  general  solu¬ 
tion  of  Eq.  (II)  in  each  layer  is  given  by 


- Z+  —Z+f* 

R  E 


cm 


where  A  and  B  are  unknown  constants.  It  is  evident  that 
there  are  four  unknowns,  two  in  each  layer.  The  B s  arc  elimi¬ 
nated  by  enforcing  the  boundary  conditions  $=0  at  z  =  ta 
and  s  =  Evidently,  we  are  assuming  zero  applied  bias 

across  the  structure.  By  pinning  the  surface  potential,  the 
surface  space  charges  Q\  and  <24  (see  Fig.  1)  brought  about 
by  the  divergence  of  the  polarization  vector  at  the  surface  are 
effectively  screened  such  that  the  electric  field  in  the  free 
space  above  and  below  the  slab  ?s  zero,  a  reasonable  and 
physically  plausible  result.  Die  relationship  between  the  two 
As  is  established  from  the  continuity  of  the  electric  displace¬ 
ment  and  is  obtained  by  integrating  Eq-  (U)  across  the  inter¬ 
face  (see  Fig.  1)  as  follows: 

C|f  +  y 

s—  =PJ|  4-+(2*3]7II+ett7:i)|  .-+^Qflau-  (13) 

^  I,-  '*  '« 

41 

From  Eqs.  (12)  and  (13)>  the  relation  between  the  As  in  the 
respective  layers  is  given  by 

AGaN=AAfcGaN+c0K2Dl  (14) 

From  Eq.  (14)  and  from  the  continuity  of  qf*  at  z-ta  =■  we  can 
solve  for  all  of  the  remaining  unknowns. 

With  the  electrostatic  potential  now  known  in  terms  of 
the*  as  yet,  undetermined  strain  tensor,  the  electric  fields  in 
the  two  layers  arc  then  given  by 


$4*  &  c 

(ID 

where  et>  is  the  electrostatic  potential  and  rz-2D  is  the  sheet 
concentration  of  an  ideal  2DEG  modeled  as  a  S  function  ai 
the  AlGaN/GaN  interface.  Equation  (11)  as  solved  analyti¬ 
cally  for  ^  subject  to  the  boundary  conditions  <£  =  0  at  z 


(J5) 


and 


^"■■-2^14+2^^-414 +4rb 


Eh.  =  -t, 


(16) 
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where  Fm=  FS(];,)~  PS{t)+^n2u  Equation  (9)  can  be 
readily  minimized  with  respect  to  y,-,  and  E.  subject  to  the 
constraints  expressed  in  Eqs,  (lO),  (15)„  and  ( 16}  using  com¬ 
mercial  symbolie-niathematkal  programs.  For  convenience, 
we  have  also  formulated  the  problem  in  the  form  of  a  matrix 
equation  of  the  form 


f  A 1 1 

A  S  l 

Au\ 

( 

^21 

A 

y'tr. 

''  ^31 

A  52 

aJ 

U J 

(17) 


in  which  the  matrix  elements  are  obtained  by  differentiating 
the  energy  functional  with  respect  to  the  unknown  strain 
components.  The  constraints  are  already  built  into  Eq.  (17) 
by  eliminating  yblX  ?  E*. ,  and  Eb.  from  Eq.  (9}.  Equation  (17) 
can  be  solved  either  symbolically  or  numerically.  We  have 
provided  the  symbolic  solution  as  an  auxiliary  item,lrt  Ap¬ 
pendix  A  lists  the  matrix  elements. 

The  fully  coupled  model  described  above  reproduces  the 
strain  tensor  from  the  standard  model  in  the  limit  that  (he 
piezoelectric  stress  tensor  i.s  set  to  zero.  The  results  are  as 
follows: 


r!* 


Frj* 


-O' 


*ba l ^  J?. ■*"  l ^3?^ 


and 


08) 


u|  „  hy(flb-aJcp 

for  the  AlGaN  layer  and 

=  _ 


and 


T%z\ 


*U* 


— -0 


>&■(  ^  b  ^  n  )  G  1 


(19) 


(20) 


(21) 


for  the  GaN  layer  where  the  constant  factor  R  is  defined  as 


( C  j  t  +  )  C  —  2  C 1 32 


(22) 


It  is  further  seen  from  Eqs.  (IS)- (21)  that  in  the  Limit 
%>tfi „  the  well-known  uncoupled  result  for  the  axial  strain  in 
an  epitaxial  AlGaN  layer  pseudomorphically  strained  on  a 
thick  GaN  buffer  is  recovered: 


r'1' 
*- 1  ■ 


7? "  I  fat  -  0  4  ^  S'* .  7*l  1 1  tt *  ' 


Ct 


33 


(23) 


In  this  limit,  the  strain  in  the  GaN  layer  is  reduced  to  zero,  at 
least  ideally,  as  the  thick  GaN  layer  serves  as  a  fixed  anchor. 
The  uncoupled  electric  held  in  the  AlGaN  layer  in  this  limit 
then  becomes 


r-ai  _  _ _  fy 


-11^31  13 1 


(24) 


We  can  also  extract  the  fully  coupled  results  for  the 
strain  and  electric  fields  in  the  limit  ti ^ta¥  i.e.  the  usual 
HFET  configuration,  and  compare  them  with  the  uncoupled 
results.  In  this  case,  fully  coupled  axial  strain  in  the  AlGaN 
layer  is  given  by 

\  c5j(*flc3j+fV>  ! 


4- 


nDCL 


(25) 


and  the  fully  coupled  electric  field  by 


r:i|  = 


C^PM 


EflC-tri  +  e 


-2 


■  31^3? 


r3SL  IS 


£-fTC“j  +  -C^ii2 


(26) 


In  Eqs.  (23)t  (24),  (25).  and  (26).  ■yj \x  takes  the  asymptotic 
limit  (ab-aa)/ua  as  lfr— The  first  term  in  Eq.  (25)  is  the 
result  from  the  standard  model.  The  other  two  terms  arc  due 
to  electromechanical  coupling. 

It  is  worth  pointing  out  that  the  foregoing  results  also 
apply  to  a  free-standing  AlGaN/GaN  superlattice  in  which 
the  slab  of  Fig.  I  forms  a  period  of  the  superlattice.  This 
outcome  follows  from  enforcing  the  pseudomorphic  bound¬ 
ary  condition  throughout  the  entire  cross  section  of  the  slab, 
i,e..  for  all  z-  In  turn,,  this  result  is  a  consequence  of  the 
simplifying  assumption  made  at  the  outset  that  ^  =  0  for 
In  addition,  the  electrostatic  potential  will  he  subject  to 
periodic  boundary  conditions  at  z=  0  and  z=tn  +  h  We  can 
assume  without  loss  of  generality  that  the  potential  in  the 
superlattice  case  can  be  set  to  zero  at  the  two  ends,  as  we 
have  done  in  the  present  case.  A  nonzero  potential  at  the  two 
ends  in  the  periodic  system  simply  represents  a  rigid  shift  of 
the  potential  distribution  function  and  does  not  change  the 
electric  field.  Thus  a  superlattice  formed  from  a  stack  of 
slabs  discussed  in  the  present  work  will  have  identical  strains 
and  electric  fields  within  each  period  as  obtained  for  our 
model  structure.  Other  free-standing  structures  for  which  the 
model  applies  include  free-standing  cantilever  microelectro- 
mechanical  systems  (MEM£)h,ia  utilizing  AIN  or  GaN.  The 
change  in  the  electrostatic  field  brought  about  by  applied 
external  stresses  on  the  MEMS  device  can  be  calculated  by 
minimizing  the  enthalpy  subject  to  the  boundary  conditions. 


III.  RESULTS  AND  DISCUSSION 

Table  I  shows  the  material  parameters1^10-14  used  in  the 
calculations.  The  signs  of  the  polarization  parameters  are  de¬ 
fined  in  relation  to  the  [0001]  direction:  a  negative  sign 
means  that  the  vector  is  in  the  [OOflf]  direction. 

For  the  model  structure  of  Fig.  Ls  we  choose 
=  300  A  and  4 =500  A.  The  300  A  AlGaN  layer  is  typical 
of  most  high-power  HFETs.  It  is  recognized  that  the  GaN 
layer  used  in  the  model  structure  is  much  thinner  than  that 
permitted  by  current  technology,  l;or  example,  using  a  laser 
lift-off  process  and  growth  via  rrvetalorganic  chemical  vapor 
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TABLE  L  Strare-rriated  material  paremewrs  used  ID  [he  present  model,  '[he 
elastic  stiffness  constants  are  in  units.  of  GPn  anti  the  piezoelectric  coeffi¬ 
cients  and  spoflhmeous  polarization  in  units  of  C/m1. 


Materia] 

C,l 

^L2 

a  (A) 

p* 

eJe0 

AIN 

3L>6J 

137“ 

16Sfl 

373h 

3.1 12* 

-0.5SJ 

1.55*1 

-  0.0S  b r 

S.51 

Crflfv 

367* 

135“ 

1 03IJ 

mx 

X\W 

-0  36' 

r 

-a.m1 

ILF 

'Reference  19. 
hRe  Terence  20, 
^Reference  2 1 . 
"Reference  22. 
cReferenas  23 . 
'Reference  15. 
Reference  24, 


deposition,  free-standing  layers  of  nitrides  ca n  be  produced 
successful ty  only  for  relatively  thick  layers25  in  the  vicinity 
of  5  /mi.  The  free-standing  layers  are  even  thicker  for  epil- 
ayers  grown  by  hydride  vapor  phase  epitaxy,  reaching  be¬ 
tween  250  and  300  jcm.26,27  The  fully  coupled  model  pre¬ 
sented  herein  is  quite  general  and,  as  shown  In  Sec.  1L  can 
readily  reproduce  the  results  for  (hick  GaN.  Figure  2  shows 
the  calculated  strain  tensor  as  a  function  of  the  A I  fraction 
using  the  fully  coupled  mode]  for  a  model  structure  that  is 
depleted,  i.e.>  «2D=0-  Following  standard  convention,  a 
positive  sign  indicates  dilation  and  a  negative  sign  contrac¬ 
tion.  Unlike  the  situation  of  a  AlGaN  layer  on  a  semi-infinite 
GaN  buffer  both  layers  are  now  strained,  with  the  in-plane 
strain  tensile  in  the  AlGaN  layer  and  compressive  in  the  GaN 
layer. 

Next  we  examine  the  deviation  between  the  fully 
coupled  and  uncoupled  models  for  a  depleted  slab.  This  de¬ 
viation  A  is  defined  as 


_  ilib)  untw|4ed 
,  a<b:i  _  iji  _ 


(27) 


where  i  is  x  or  z*  the  superscript  ‘V*  or  "b”J  refers  to  a 
particular  layer,  and  yfi  is  given  by  the  solution  to  Bq.  (17) 
for  the  coupled  case  and  by  Eqs.  ( 1 8)-{2 1 )  for  the  uncoupled 
case.  The  calculated  results  are  shown  in  Figs,  3  and  4  for  the 
in- plane  and  out-of-plane  strains,  respectively.  It  is  seen  from 


FSC,  3.  Eteviaiion  of  the  uncoupled  from  the  fully  coupled  in-plaiw  strain  in 
the  bit&yer  Stab  with  n2[>=G, 

Fig.  3  that  the  deviation  for  is  quite  small.  A  possible 
reason  is  that  the  electromechanical  coupling  of  yxx  into  the 
equation  of  state  [Bq.  (7)]  occurs  through  the  piezoelectric 
coefficient  e3!  which  is  much  smaller  than  e3i,  as  seen  from 
Table  1.  The  electromechanical  coupling  of  y_..  into  the  equa¬ 
tion  of  state  is  through  e 33  which  is  quite  large,  especially  for 
AIN.  Jt  is  dear  that  the  error  in  the  calculated  strain  along  the 
growth  axis  is  quite  significant,  reaching  30%  in  the  AlGaN 
layer  for  a  mole  fraction  of  0.3.  This  error  will,  in  turn,  have 
a  significant  impact  on  the  calculation  of  electronic  and  op¬ 
tical  properties  that  depend  on  accurate  knowledge  of  the 
strain.  Some  examples  of  such  quantities  include  the  eigen¬ 
states  and  the  complex  dielectric  function  of  the  slab. 

So  far  it  has  been  assumed  that  the  slab  is  depleted  with 
the  only  charge  present  being  a  fixed  space  charge  P  given 
by 

p—  p£iiii_  pp{ b>_  ppi* j  j  n g ) 

where  P^  and  PF^  are  the  piezoelectric  polarizations 

<»> 

and 

F*b’=  2e^+eb^.  (30) 


FIG.  2.  Calculated  strain  for  the  SKuaare  of  Fig.  3  as  a  function  of  Ihc  At 
fraction  with  f,  =  3D(jAb  f*=500A,  and  n2D=0,  using  the  hilly  coupled 
model.  The  four  sirain  components  are  denoted  in  Fig.  L 


FIG,  4.  Deriaiion  of  the  uncoupled  from  the  fully  coupled  our-of-planc 
strain  in  the  bitayer  slab  with  n1D=0. 
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FIG.  5  Calculated  ouhof-planc  strain  in  Lhe  AlGaN  layer  for  (he  standard 
and  coupled  models.  TWO  coupled  cases  arc  shown,  one  without  frvt  elec¬ 
trons  and  one  wilh  free  elect™*;. 


For  a  cation-faoed  structure,  the  standard  growth  orientation 
for  most  HFBTs,  P  is  positive;,2-'  where  P  =  Q2~  Qy  in  Fig. 
1.  Depending  on  the  growth  and  surface  conditions,  donors 
may  be  present  and  may  contribute  free  electrons  to  the  in¬ 
terlace,  Under  certain  conditions,  the  space  charge  may  be 
partially  screened  by  a  2PEG-  For  instance,  using  our  own 
self-consistent  Schrodinger- Poisson  model,"1  we  find  that 
the  2 DEG  is  almost  90%  of  the  magnitude  of  the  fixed  space 
charge  for =  300  A  and  typical  HFET  AS  mole  fractions, 

Tbe  exact  origin  of  the  2DEG  is  not  well  known.  From 
basic  electrostatics*  the  sum  of  the  polarization-induced 
space  charges  must  be  zero.  This  means  that  Qi  +  Qi  + 

+  Qi  =  0  (see  Fig.  I).  There  remains  die  issue  of  how  to 
account  for  the  2DEG  widiin  a  theoretical  framework  while 
ensuring  that  the  charge  neutrality  law  is  satisfied.  In  our 
numerical  charge-control  model2*  we  assume  that  the  2 DEG 
originates  from  surface  donors.  Charge  balancing  is  then  ef¬ 
fected  by  the  2DEG  being  equal  to  the  sheet  density  of  ion¬ 
ized  donors.  Other  possible  sources  of  the  2DEG  include 
unintentional  dopants  and  deep-level  traps.  The  magnitude  of 
the  2DEG  will  depend  on  [he  band  lineup,  the  thickness  of 
the  layers,  and  the  Al  fraction  of  the  barriers.  In  addition, 
there  is  both  a  minimum  thickness  tn  and  a  minimum  Al 
fraction  for  the  2D  EG  to  form,  as  showrn  in  a  recent  numeri¬ 
cal  work.’0  In  the  present  wrork,  however,  the  formation 
threshold  of  the  2  DEG  as  a.  function  of  the  A  l  fraction  is  not 
taken  into  account.  It  has  also  been  shown  numerically^0  that 
the  2D£G  density  is  a  linear  function  of  the  A!  fraction  for  a 
given  ia  .  We  make  the  latter  assumption  in  the  present  work. 
Further  we  assume,  conservatively,  that  nS[J=O.EP/^0  and 
examine  the  effect  of  free -electron  screening  on  the  electro¬ 
mechanical  coupling. 

Figures  5  and  6  show-  ytz  and  y!?.  „  respectively,  calcu¬ 
lated  using  the  standard  and  coupled  models.  The  in-platne 
strains  are  not  shown,  since  the  effect  of  electromechanical 
coupling  on  the  in -plane  strain  components  is  weak.  It  is 
evident  from  the  calculated  strains  that  the  effect  of  coupling 
is  to  reduce  the  magnitude  of  the  strains  relative  to  those  of 
the  standard  model,  especially  when  the  structure  is  depleted. 
This  result  is  not  surprising,  since  we  would  expect  the  pi¬ 
ezoelectric  fields  to  induce  forces  to  oppose  any  applied 


FIG.  6.  Calculated  out-of-plaue  strain  in  the  GaN  layer  for  the  standard  and 
coupled  models.  Two  coupled  cases  are  shown,  MIC  without  free  electrons 
and  one  with  free  electrons. 


forces  such  as  those  due.  to  the  pseudomorphic  interface  con¬ 
dition.  Another  result  that  is  evident  from  Figs.  5  and  6  is 
that  the  2DEG  screening  brings  the  strains  of  the  coupled 
model  closer  lo  those  of  the  standard  model.  As  will  be  seen 
shortly,  the  effect  of  screening  is  to  reduce  the  built-in  elec¬ 
tric  holds.  This,  in  turn,  will  reduce  the  electromechanical 
coupling,  as  seen  from  Eq,  (7). 

Figures  7  and  &  show-  the  calculated  electric  fields  in  the 
AlGaN  and  GaN  layers,  respectively,  using  the  standard  and 
coupled  models-  These  fields  help  to  explain,  in  part,  the 
large  deviation  of  the  standard  strain  fields  from  their 
coupled  counterparts.  As  an  example,  for  a  Al  fraction  of  0.3, 
the  electric  held  in  the  AlGaN  layer  ls  about  2  MV/cm*  pro¬ 
vided  that  the  structure  is  depleted.  Such  a  large  held,  if 
present,  will  give  rise  to  a  strong  electromechanical  coupling 
and,  therefore,  a  largo  discrepancy  between  the  standard  and 
coupled  models.  When  there  is  a  2 DEG,,  however,  The  high 
space  charge  density  that  gives  rise  to  the  electric  held  is 
partially  neutralized,  depending  on  the  magnitude  of  the 
2 DEG.  For  fr2D=Q.£P,  the  field  in  the  AlGaN  layer  is  re¬ 
duced  to  about  Q.4  MV/cm  for  a  A l  fraction  of  03-  Evidence 
for  the  smaller  field  is  seen  in  recent  phntoreflectanee 


FIG.  7.  Calcidaled  electric  field  in  the  AlGaN  layer  fur  the  standard  and 
fully  coupled  models  showing  rite  effect  of  free-electron  screening  in  each 
cast. 
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FIG.  E.  Calculated  electric  field  in  the  GaN  Layer  for  the  standard  and  Fully 
coupled  models  showing  the  effect  of  finee-electron  screening  in  each  case. 

measure ments'3  -2  on  AlGaN/GaN  beterostructures  in  which 
barrier  fields  in  the  vicinity  of  0.3  MV/cm  were  reported  for 
250  A  thick  barriers  with  A!  fractions  in  rhe  range  0.06-0.1. 
An  electric  field  of  opposite  sign  occurs  in  the  GaN  layer 
with  similar  trends  between  the  coupled  and  standard  mod¬ 
els.  Its  magnitude  is  decreased  with  the  thickness  of  the  GaN 
layer.  It  is  recognized  that  the  calculation  of  the  electric  field 
in  the  GaN  layer  is  more  involved  than  the  classical  model 
presented  allows  here.  The  formation  of  the  notch  in  the 
conduction  hand  edge  in  which  the  2DEG  resides  is  repro¬ 
ducible  only  wdthin  the  scope  of  a  quantum -mechanical  cal¬ 
culation.  In  spite  of  approximations  used,  the  calculated 
fields  shown  here  appear  to  be  quite  plausible. 

IV,  SUMMARY  AMD  CONCLUSIONS 

Tn  summary,  a  fully  coupled  electromechanical  model 
has  been  presented  for  the  strain  and  electric  fields  in  a  free¬ 
standing  hi  layer  AiGaN/GaN  slab.  The  electric  enthalpy  of 
the  stab,  with  contributions  from  the  piezoelectric  and  spon¬ 
taneous  polarizations,  is  minimized  subject  to  [he  constraints 
imposed  by  the  pseudomorphic  interface  condition  and  the 
Poisson  equation.  Closed-form  expressions  for  the  strain 
components  and  electric  field  in  each  layer  are  obtained.  The 
results  for  the  bilayer  system  are  also  appropriate  for  free¬ 
standing  binary  superlattices  with  negligible  bowing.  Large 
discrepancies  between  the  standard  and  coupled  models  are 
found  tor  depleted  structures  due  to  large  built-in  electric 
fields.  The  electric  fields  are  reduced  when  the  polarization- 
induced  space  charge  is  screened  by  a  2DEG  As  a  conse¬ 
quence,  the  discrepancy  between  the  standard  and  coupled 
models  is  reduced  when  a  2DEG  is  present 
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APPENDIX:  MATRIX  SOLUTION  OF  THE  STRAIN 
FIELD 

The  matrix  elements  required  to  solve  Eq.  (17)  are  given 
in  this  Appendix.  As  before,  ‘V  refers  to  the  ASGaN  layer 
and  “b1"  to  the  GaN  layer 

A  I  j  G^fjl  ^1  I  ""f  ^  l2.}(  ^-fr^fr)  f  d (V ra1® b  1 

XCC^,  +  cbn)  +  2  tM*  31a«  -'3lfli)a-  (Al) 

(A2) 

a  ] J=  tba„ai,C'lsUaBb+  tbsa)  - <3^  J- 

(A3) 

A  7,  flCI h  "t"  „)  4  2/^1  ^33(^3  \Ctfr  £ 

(A4) 


^22” *fr£n)  tn^bat?€^  * 

(A5) 

A  tj  =  —  iai  &&  i>e  3  jt,  a 

(A6) 

A  ji  =  -  e  +  2tlj)„Cb]  J, 

(A7) 

A  JJ—  —  t„I  iflfrS  3  J , 

(AS) 

B I  =  tatb(e  Jiflfr-  e  FMb)  -  PSi*'  +  ^id) 

(A9) 

^(taeb+tbe«h 

(A10) 

X(PSibi-PSw+eQnllt), 

(All) 

and 

Bi=~  tJba^P^-  PSw+  e<,n2a)  +  2{ab -  a  J 

(A12j 
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